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Abstract 

The Two-dimensional Implicit Thermal Response and Ablation program, TITAN, 
was developed and integrated with a Navier-Stokes solver, GIANTS, for multi- 
dimensional ablation and shape change simulation of thermal protection systems in 
hypersonic flow environments. The governing equations in both codes are discretized 
using the same finite-volume approximation with a general body-fitted coordinate 
system. Time-dependent solutions are achieved by an implicit time marching technique 
using Gauess-Siedel line relaxation with alternating sweeps. 

As the first part of a code validation study, this paper compares TITAN-GIANTS 
predictions with thermal response and recession data obtained from arc-jet tests recently 
conducted in the Interaction Heating Facility (IHF) at NASA Ames Research Center. The 
test models are graphite sphere-cones. Graphite was selected as a test material to 
minimize the uncertainties from material properties. Recession and thermal response data 
were obtained from two separate arc-jet test series. The first series was at a heat flux 
where graphite ablation is mainly due to sub imation, and the second series was at a 
relatively low heat flux where recession is the result of diffusion-controlled oxidation. 
Ablation and thermal response solutions for both sets of conditions, as calculated by 
TITAN-GIANTS, are presented and discussed in detail. Predicted shape change and 
temperature histories generally agree well with :he data obtained from the arc-jet tests. 

Nomenclature 

A = area, m 2 

= ml p e u e C m , dimensionless mass blowing rate 

B a = pre-exponential constant in Eq.(6), s' 1 
Ch = Stanton number for heat transfer 
C m = Stanton number for mass transfer 
c p = specific heat, J/kg-K 
E = total energy per unit volume, J/m 3 
E a = activation energy in Eq.(6), J/kmol 
F = radiation view factor in Eq.(l 1) 
h = enthalpy, J/kg 

h = enthalpy of the partial-pyrolyzed solid, defined in Eq.(4), J/kg 
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hj W = enthalpy of species j respect to surface temperature in Eq.(l 1), J/kg 

H r = recovery enthalpy, J/kg 

k = thermal conductivity, W/m-K 

m = mass flux, kg/m -s 
2 

p = pressure, N/m 

qc = conductive heat flux, W/m 2 

qR = radiative heat flux, W/m 2 

q v = heat flux due to species diffusion, W/m 2 

R = universal gas constant, J/kmol-K 

R n = nose radius, m 

T = temperature, K 

t = time, s 

u = fluid velocity, m/s 

v = local grid velocity, m/s 

x = Cartesian coordinate system, m 

Z* = coefficient in Eq.(l 1), defined in Ref. 14 

a = surface absorptance 

e = surface emissivity 

r = volume fraction of resin 

X = blowing reduction parameter in Eq.(12) 

p = total density, kg/m 3 

p 0 = original density, kg/m 3 

p r = residual density, kg/m 3 

a = Stefan-Boltzmann constant, W/m 2 -K 4 

T = mass fraction of virgin material defined in Eq.(3), or shear stress in Eq.(9) 

¥ = decomposition reaction order in Eq.(6) 

£,r|= general body fitted coordinate system defined in Fig.2 

v = species diffusion flux, m/s 

w = species source term in Eq.(8), kg/m 3 -s 

V = gradient, m' 1 

subscripts 

c = char 

e = boundary-layer edge 
g = pyrolysis gas 

i = material density component (A, B, and C), or direction component 
j = surface species, or direction component 
s = gas species or stagnation point 
v = virgin 
w = wall 

Introduction 

Heatshields of spacecraft for planetary mission typically use thermal protection 
system (TPS) materials which pyrolyze and ablate at high temperature for mass-efficient 
rejection of aerothermal 
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heat load. Pyrolysis is an internal decomposition of the solid which releases gaseous 
species, whereas ablation is a combination of processes which consume heatshield 
surface material. A reliable numerical procedure that can compute surface recession rate, 
mass loss, in-depth pyrolysis, and internal temperature time histories under general 
heating conditions is essential for the design and sizing of ablating spacecraft TPS 
materials. 

The purpose of this paper is to compare TITAN-GIANTS predictions with arc-jet 
data and to study the accuracy of the simulation. Thermal response and recession data 
were obtained from arc -jet tests recently conducted in the Interaction Heating Facility 
(IHF) at NASA Ames Research Center. Graphite ablation and thermal properties have 
been widely studied, so graphite was selected as the test material to minimize the 
uncertainties from material properties (the accuracy of shape change and temperature 
history predictions can thus be better determined). However, the modeling of pyrolysis 
gas cannot be examined from graphite data. The pyrolysis gas effect will be studied in 
future work. 

The recession and thermal response data studied here are from two separate arc- 
jet test series. The first series was at relatively high heat flux in which graphite ablation 
was mainly due to sublimation, and the second series was at relatively low heat flux in 
which the recession was the result of diffusion-controlled oxidation. Typical B' vs. 
temperature curves for graphite in air at various surface pressures (0.01, 0.1, and 1 atm) 
are shown in Fig. 1. Above 3000 K, B' is a strong function of temperature, and 
sublimation is the primary ablation process. At lower temperature, B' is relatively 
constant (approximately 0.175), and diffusion-controlled oxidation determines the 
ablation rate. At temperatures below 1 100 K, the ablation is kinetically controlled and is 
beyond the scope of this work. 

Background 

The Charring Material Thermal Response and Ablation Program 1 (CMA) was 
developed by Aerotherm Corporation in the 1960s. This technique solved the one- 
dimensional internal energy balance and decomposition equations coupled with the 
ablating surface energy balance condition to simulate the response of ablative heatshields 
in hypersonic flows. The Fully Implicit Ablation and Thermal Response program 2 (FIAT) 
was developed at NASA Ames Research Centei in the 1990s. Both FIAT and CMA solve 
the same governing equations. FIAT, which is numerically more stable, solves a wider 
range of problems compared with CMA, and has been used for TPS sizing calculations 
for various NASA space missions, such as Mars Pathfinder 3 and Stardust. 4 

The Two-dimensional Implicit Thermal Response and Ablation (TITAN) 
program 5 was developed and coupled with the GIANTS (Gauss-Seidel implicit 
aerothermodynamic Navier-Stokes code with thermochemical surface boundary 
conditions) code 6 ' 9 using a loosely coupled method to perform high fidelity thermal 
response and shape change simulation for charring materials. The governing equations in 
both TITAN and GIANTS are discretized usirg the same finite-volume approximation 
with general body fitted coordinates. The solution is achieved by an implicit time 
marching technique using Gauess-Siedel line relaxation with alternating sweeps. A fluid- 
solid coupled simulation is required to accurately predict shape change and aerothermal 
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heating, because aerothermal heating distribution is very sensitive to body shape, and 
shape change is a strong function of surface heating distribution. TITAN-GIANTS was 
also integrated with the commercial finite-element code, MARC, 10 to perform thermal 
and structural analysis on the complex geometries of actual vehicles. Details of 
TITAN’s development and its formulation hr.ve been discussed in a previous paper. 
Application of TITAN-GIANTS and TITAN-MARC for full body heat-shield/structure 
analysis of a Mars sample return Earth Entry V ehicle and for simulation of a flat-faced 
cylinder arc -jet test model are demonstrated in Reference 1 1 . 

Arc- Jet Facilities 

Arc -jet facilities are used to simulate the aerothermal heating environments 
experienced by space vehicles during atmospheric entry. The arc -jet tests discussed in 
this paper were conducted in the Ames 60 MW Interaction Heating Facility. In the IHF, 
air is heated by the electric discharge in an 8-cm diameter, water-cooled constricted arc 
column. The hot gas is then expanded through a converging-diverging nozzle into an 
evacuated test chamber where the test model is located. The heat flux and pressure on the 
model surface can be varied by changing the el metrical power input and the configuration 
of the nozzle, as well as the model location from the nozzle exit. The test model is 
mounted onto a mechanical swing arm and positioned between four and fourteen inches 
from the nozzle exit. The swing arm inserts the test model into the gas stream when the 
desired test conditions are reached, and then retracts the model from the stream after the 
desired exposure time has elapsed. 

Two-dimensional Implicit Thermal and Ablation Program - TITAN 

The ablating heat-shield material thermal response and shape change 
computations are performed using TITAN. The governing equations include energy 
conservation and a three-component decomposition model. The surface energy balance 
condition is solved with a moving grid to calculate the shape change due to surface 
recession. The internal energy balance is a transient thermal conduction equation with 
additional pyrolysis terms 5 

p Cp — = V ■ (k VT) - ( h i g - /T)V • m s + m K • V h g + vp c P VT (1) 

d t 

The individual terms in Eq. 1 may be interpreled as follows: rate of storage of sensible 
energy, net rate of thermal conductive heat flux, pyrolysis energy-consumption rate, net 
rate of energy convected by pyrolysis, and convection rate of sensible energy due to 
coordinate system movement. The local specific heat is formulated from functions of 
temperature input for both virgin material and char. In partially pyrolyzed zones (p c < p < 
p v ), the specific heat is obtained from the mixing rule 


c P = r c pv +(l-T)c pc , V) 

where t is the mass fraction of virgin material in a mixture of virgin material and char: 
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( 3 ) 


T _ i-Pc/P 

1-Pc/Pv ' 

The thermal conductivity, k, is weighted in the same manner. The pyrolysis gas 
enthalpy, h g , is an input function of temperature and pressure. The enthalpy of the partial- 
pyrolyzed solid, h , is defined as 

h _ PA~PA (4) 

Pv-Pc 


A three-component decomposition model is used. The resin filler is assumed to 
consist of two components which decompose separately, while the reinforcing material is 
the third component which can decompose. The instantaneous density of the composite is 
given by: 

p=rxp, + p,)+(i-r)p,,, (5) 


where A and B represent components of the resin, and C represents the reinforcing 
material. T is the volume fraction of resin and is an input quantity. For graphite, the 
volume fraction, T, is 0, and the reinforcing material density, p r , is the graphite density. 
Each of the three components can decompose by the relation: 


5P, 

dt 


. — Eai . 

= -B at exp(-^ r )p 0 , ( 


— — ——) + v Vp , 


( 6 ) 


where p ri is the residual or terminal density of component i , and p m is the original 

density of component / ( / = A, B, and C). The' motion of pyrolysis gas is assumed to be 
quasi-steady and one-dimensional (along the i] direction), and thus the mass flow rate of 
pyrolysis gas at the surface is calculated using the following approximation: 



(7) 


The governing equations are discretized using a finite-volume approximation with 
a general body fitted coordinate system. The time-dependent solution is achieved by a 
implicit time marching technique using Gauess-Siedel line relaxation with alternating 
sweeps. If time accuracy is not of interest in the computation, the time step should be as 
large as possible, and the number of alternating sweeps in each time step should be 1 . A 
time-accurate solution can be achieved by increasing the number of alternating sweeps in 
each time step or by reducing the time step. The: computational grid is compressed during 
the course of computation to reflect the surface recession. Further details of the numerical 
procedures are presented in Reference 5. 
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Navier-Stokes equation solver - GIANTS 


To estimate the aerothermal heating distribution over a large angle blunt body, the 
hypersonic flow simulation is performed using the Navier-Stokes solver, GIANTS. The 
GIANTS code solves the time-dependent conservation equations of mass, momentum, 
and energy for the chemical and thermal non- equilibrium flow-field. The species mass 
conservation equation is given by 12 


dt dXj 


(P<Uj) = --^-(PsV sj ) + w,, 


( 8 ) 


the total momentum conservation is written as 


_d 

dt 


(P«,) + 



dr y 

(pUiUi> 


(9) 


and the total energy conservation as 

dE d „ 3 , 9 , 

— + — ((£ + p)Uj ) = -— (qj +q VJ )-^-(u,x, J , 


The governing equations are discretized using the finite-volume method. The numerical 
method (fully implicit and Gauss-Seidel line relaxation) used to solve the discretized 
equations is exactly the same as that used in the TITAN code. This technique has been 
shown to yield steady-state results efficiently. A five species (N 2 , O 2 , NO, 0,and N) gas 
model is used to accurately model the high temperature behavior of air. A bifurcation 
diffusion model 13 was also implemented in this version of the code to correctly compute 


multi-species mass diffusion. 

TITAN and GIANTS are interfaced through the surface energy balance boundary 
condition. The conditions at the ablating surface are determined by convective and 


radiative heating and by surface thermochemical interactions with the boundary layer 
gases. The surface energy balance equation employed is of the convective transfer 
coefficient type. This energy balance equation takes the following form 


p e u c C h (Hr-h e ») + P'U'Cj'ZiZ], -Z] w )h'; - B' hw\ + rhchc + m g h g 

+ a w q Rv -FoeT*-q CM =0 


( 11 ) 


The first term in Eq. 1 1 represents the sensible convective heat flux. The sum of the 
second, third, and fourth terms in Eq. 1 1 is defined as the total chemical energy at the 
surface. The Z* terms represent transport of chemical energy associated with chemical 
reactions at the wall and in the boundary layer. 1 ' 1 The Z* driving forces for diffusive mass 
transfer include the effects of unequal diffusion coefficients. The fifth and sixth terms are 
the radiative heat fluxes absorbed and re-radial ed by the wall, respectively, and the last 
term, qc w , represents the rate of conduction into the material. The radiative heat flux, 
q Rw , such as radiation from a hot shock layer, is an input quantity. The prediction of 

shock layer radiation is not currently included in the coupled TITAN-GIANTS 
simulation. Here g' is the normalized mass blowing rate. The tables of g' for charring 
materials can be pre-generated using ACE 14 or MAT 15 assuming chemical equilibrium. 

The coupled TITAN-GIANTS simulation does not include surface mass blowing and 
ablation chemical species in the Navier-Stokes computation. Instead, a blowing 
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correction parameter is introduced to account for the reduction in transfer coefficients due 
to the transpiration effect of the mass injection into the boundary layer. This approach 
was taken because previous studies of heat-shield design for Mars Pathfinder and 
Stardust 4 , demonstrated that the computational time increased by about at least one order 
of magnitude if a blowing surface condition and ablation chemical species with their 
associated reactions were included in a flow simulation. In addition, the reaction rates of 
many gas phrase chemical reactions are unknown or not well studied. Thus, the 
prediction of Navier-Stokes solver becomes unreliable because of the uncertainties on the 
gas phrase chemical reaction rates. The introduction of a blowing reduction parameter 
can give a reasonable prediction for many applications and significantly cut the CPU 
time. The blowing rate correction for convective heat transfer is 

C h _ ln(l + 2A5') 

2 XB' 

where A is the blowing reduction parameter, C h is the heat transfer coefficient for the 
ablating surface, and C M is the heat transfer coefficient for a nonablating surface. When 

A = 1 / 2 , Eq.( 1 1 ) reduces to the classical laminar-flow blowing correction. 1 6 

The computation starts with the TITAN code using initial estimates of surface 
heat flux and pressure without shape change as the surface boundary condition. At each 
time increment, TITAN obtains the in-depth thermal response solution based on current 
boundary conditions, and also calls the flow environment code to update the front-face 
energy balance boundary condition if necessaiy. The non-ablating surface heating and 
pressure are calculated by GIANTS, and a blowing reduction parameter of 0.5 is used in 
TITAN to take into account the laminar flow blockage due to surface mass blowing. 
When the maximum local surface recession exceeds a predefined criterion since the last 
surface convective heating was updated, a flow-field grid is generated based on the 
current body geometry, and then the flow simulation routine is called upon to compute 
the new aerothermal heating environment. Each call to the flow environment routine is a 
steady state calculation. As expected, the GIANTS calculation is much more 
computationally intensive than the thermal response computation performed by TITAN. 
For a TITAN-GIANTS coupled simulation, most of the CPU time is consumed by the 
flow environment calculation. Thus, the CPU tune required for a simulation is primarily 
determined by the efficiency of the GIANTS code. 

Results and Discussion 

The first series of arc -jet tests were conducted with a stream enthalpy of 27 
MJ/Kg, a stagnation point cold-wall heat flux of 2100 W/cm 2 , and a stagnation point 
pressure of 0.75 atm. The length of the heat pulse is 30 sec. The Poco graphite model is a 
10° half angle sphere-cone with nose radius of 1.905 cm (0.75”). The thermal properties 
of Poco graphite are from Ref. 17. The total length of the model is 8.89 cm (3.5”). The 
computational grids for both solid and fluid are shown in Fig. 2, where the solid grid is 
divided into two zones, and the fluid grid has a single zone. The first zone in the solid 
will be reconstructed as the recession exceeds preset criteria, and the second zone 
remains unchanged through the entire computation. Surface heat flux and pressure 
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distributions are recalculated using GIANTS when free stream conditions or the first zone 
in the solid grid is changed. The thermal response computation of the solid is transient, 
and the flow field is assumed to reach steady state instantaneously. 

The heat flux distributions over the solid surface at t = 0 and 30 sec are shown in 
Fig. 3. The stagnation point heating rate at the end of heat pulse (t = 30 sec) is about 90% 
of that at the beginning (t = 0 sec). However, the heating at 45° from stagnation point 
becomes slightly higher toward the end of the heat pulse. Surface recession tends to blunt 
the spherical nose, and thus reduce stagnation point heating. Because of the uncertainties 
in estimated arc-jet free stream conditions and material surface conditions, the GIANTS 
solution at t = 0 does not exactly match the measured stagnation-point cold- wall heat flux 
(2100 W/cm 2 ). To account for these uncertainties, the normalized heat flux distributions 
from GIANTS were scaled to exactly match the measured initial stagnation point heat 
flux. 

A comparison of predicted and measure d surface recession is shown in Fig. 4. The 
symbols for the graphite surface recession data are measured every 5 seconds at the 
stagnation point and at 45° from the stagnation point. (Details of the measurements are 
described in Ref. 18.) The dark solid lines are the predicted recession using a coupled 
TITAN-GIANTS simulation. The total recession for this 30-sec heat pulse at stagnation 
point is about 0.33 cm and at 45° from stagnation point is about 0.165 cm. The 
computation slightly over-predicts the recession at the stagnation point and slightly 
under-predicts the recession at 45° . The max mum difference between prediction and 
measurement is less than 1 0 %. 

For comparison purposes, an additional FIAT computation was performed by 
assuming that the cold- wall heat flux was 2100 W/cm 2 and the thermal diffusion was 
one-dimensional. FIAT predicts much higher recession compared with TITAN at the 
stagnation point, as shown by the lighter solid line in Fig. 4. This result is not surprising, 
because the FIAT calculation does not take into account multidimensional heat 
conduction and shape change effects. 

The B' histories at the stagnation point and at 45° are presented in Fig. 5. In the 
first 2 seconds, the stagnation point is in the diffusion-controlled regime, but as the 
surface temperature increases the ablation enters the sublimation regime, and B' levels 
off near 0.4 until the end of heat pulse. At 45° from the stagnation point, the surface 
remains in the diffusion controlled regime unti about 5 seconds, and then B' gradually 
increases to about 0.3 at 30 sec. 

The predicted surface temperature his :ory at 45° from the stagnation point is 
shown in Fig. 6 and compared with temperature measurements from an infrared (IR) 
camera and a pyrometer. At t = 30 sec, the predicted surface temperature is 3180°C , the 
temperature measured from IR camera is 3030 °C , and that from the pyrometer is 2940 
°C . The difference between the prediction and the IR is about 5% and between the 
prediction and the pyrometer is about 8 %. Further study is required to fully understand 
why the calculated temperature is higher than the measured data, and why the 
temperature reading from the IR camera is higher than that from the pyrometer. 

There were no stagnation point surface temperature and in-depth temperature 
measurements in this test series. The predicted temperature contours of the graphite 
model at t = 30 sec are presented in Fig. 7. The stagnation point temperature is above 
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3500 K, and the back face of test model, which is connected to a water cooled holder, is 
assumed to remain at 312 K. The dark solid line represents the initial surface contour. 
The area between the dark solid line and color contours is the material ablated during the 
30 sec heat pulse. 

Figure 8 presents the flow field predictions at 0 and 30 sec. The upper half is the 
flow temperature contours at the beginning of test (t = 0 sec), and the lower half is that at 
the end of the heat pulse (t = 30 sec). The Knudsen number is about 0.023. As expected, 
there is no significant change in flow pattern. The maximum fluid temperature around the 
shock is above 10,000 K. The shock stand-off distance at t = 30 sec is slightly larger than 
at t = 0 sec. This is because ablation blunts the spherical section, and consequently 
reduces the heating in the stagnation region. 

The second arc -jet test series was conducted with a stream enthalpy of 19 MJ/Kg, 
a stagnation point cold wall heat flux of 593 W/cm , and a stagnation point pressure of 
0.05 atm. The duration of heat pulse was 60 se:. The ATJ graphite composite model is a 
15° half angle sphere-cone with nose radius of 0.95 cm (0.375”). The total length of 
model is 10.16 cm (4”). The material map and computational grids for both solid and 
fluid are shown in Fig. 9. The solid grid is again divided into two zones, and the fluid 
grid has a single zone. The first 3.81 cm (1.5”) of the model is a mushroom-shaped piece 
of ATJ graphite. 19 The rest of the model is built from a AETB 35/12 insulator with a 
small region of TUFROC to prevent overheating of the bond-line to the AETB. 

The heat flux distributions over the solid surface at t = 0 and 60 sec are shown in 
Fig. 1 0. The stagnation point heating rate at the end of heat pulse (t = 60 sec) is about 
95% of that at the beginning (t = 0 sec), and again (similar to the first series) at 45° from 
stagnation point the heat flux becomes slightly higher at the end of heat pulse. However, 
the difference in surface heating between the beginning and the end of heat pulse is 
relatively small, compared with the first test series, because the total recession is lower in 
this second series of tests. 

The stagnation point recession, measured at t = 25 and 60 sec, is compared with 
predictions in Fig. 11. The solid line is the prediction from TITAN, and the dashed line is 
from FIAT (a one-dimensional calculation with constant heat flux of 593 W/cm 2 ). As 
expected, FIAT predicts slightly higher recession compared with TITAN. The maximum 
difference between TITAN and FIAT is about 9% at t = 60 sec. The measured stagnation 
point recessions at 25 and 60 sec are 0.068 arid 0.183 cm, respectively. The predicted 
recession from TITAN at 25 sec and 60 sec is 0.071 cm and 0.174 cm, respectively. The 
maximum difference between the TITAN prediction and the measurement is about 5%. 
In this test series, graphite ablation is mainly due to diffusion-controlled oxidation. Thus, 
B' remains relatively constant (approximately 0.175) through the entire heat pulse. The 
total mass loss in this second test series is smaller than in the first one. 

Three thermocouples (designated TCI, TC2, and TC3) are located on the model 
axis at depths of approximately 3.18, 4.13, and 6.67 cm from the initial stagnation point. 
The predicted temperature histories at these three locations from TITAN (dark solid 
lines) are presented and compared with the TC data in Fig. 12. Generally speaking, the 
TITAN predictions are in good agreement with .he data. The TITAN predicted maximum 
temperatures of TCI and TC2 are slightly below the data, but the predicted cooling rates 
for these two TC’s are slightly lower than data. At the location of TC3, the TITAN 
prediction is slightly lower than the data most of the time. 
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Surface temperature data are not available for this test series. Predicted 
temperature contours at t = 60 sec (end of the heat pulse) and at 200 sec (during 
cooldown) are shown in Figs. 13 and 14, respectively. The peak temperature at 60 sec is 
above 2700 K around stagnation point, and at 200 sec the hottest location is about 1200 K 
near the interface between the graphite and the AETB. The back face of model is 
connected to a water cooled holder and assumed to be at 312 K. The dark solid line again 
represents the initial surface contour. 

Figure 15 presents the flowfield predict ons at 0 and 60 sec. The upper half is the 
temperature contours at the beginning of heat pulse (t = 0 sec), and the lower half is that 
at the end of heat pulse (t = 60 sec). The Knudsen number for this case is about 0.1. 
These results are similar to the first series (Fig. 8) in that there is no significant change in 
flow pattern. The maximum fluid temperature around the shock is above 7000 K. The 
shock stand-off distance at t = 60 sec is almost the same as that at t = 0 sec. 

Conclusions 

Predictions from the integrated TITAN-GIANTS code were compared with 
thermal response and recession data obtained from arc -jet tests conducted in the 
Interaction Heating Facility at NASA Ames Research Center. The test models are 
graphite sphere-cones. The arc-jet data were taken from two separate test series. The first 
test series (q s = 2,100 W/cm 2 , and P w = 0.75 atm.) had graphite ablation in the 
sublimation regime, and the second test series (q s = 593 W/cm 2 , and P w = 0.05 atm) had 
ablation in the diffusion controlled oxidation regime. Ablation and thermal response 
solutions calculated by TITAN-GIANTS for both sets of test conditions agree reasonably 
well with the measurements. The calculations slightly over-predict the surface 
temperature at 45° from the stagnation point for the sublimation case, and slightly under- 
predict in-depth maximum temperatures for ciffusion-controlled case. The maximum 
difference in recession between TITAN predict ions and measurement is less than 1 0 %. 
A one-dimensional FIAT computation can be used to roughly estimate the stagnation 
point recession, but not to predict in-depth thermal response. 
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Figure 1 : Graphite B' curves in ai r (generated using JANAF data) 


q s = 2,100W/cm 2 
P = 0.75 atm 



Figure 2: Computational grid for test series 1 
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Figure 7: Temperature contours for graphite model at t = 30 sec 
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Figure 8: Flowfield temperature contours at t = 0 and 30 sec. 
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q s = 593 W/cm 
H = 1 9 MJ/Kg 
P = 0.05 atm 








Figure 12: In-depth temperature histories along axis at three depths 
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Figure 13: Temperature contours for composite model at t - 60 sec 



Figure 14: Temperature contours for composite model at t - 200 sec 
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